Quantum discord of two-qutrit system under quantum-jump-based feedback control
Wang Chang, Fang Mao-Fa
Synergetic Innovation Center for Quantum Effects and Application, Key Laboratory of Low-dimensional Quantum Structures and Quantum Control of Ministry of Education, School of Physics and Electronics, Hunan Normal University, Changsha 410081, China

 

† Corresponding author. E-mail: mffang@hunnu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant No. 11374096).

Abstract

This paper studies quantum discord of two qutrits coupled to their own environments independently and coupled to the same environment simultaneously under quantum-jump-based feedback control. Our results show that spontaneous emission, quantum feedback parameters, classical driving, initial state, and detection efficiency all affect the evolution of quantum discord in a two-qutrit system. We find that under the condition of designing proper quantum-jump-based feedback parameters, quantum discord can be protected and prepared. In the case where two qutrits are independently coupled to their own environments, classical driving, spontaneous emission, and low detection efficiency have negative effect on the protection of quantum discord. For different initial states, it is found that the evolution of quantum discord under the control of appropriate parameters is similar. In the case where two qutrits are simultaneously coupled to the same environment, the classical driving plays a positive role in the generation of quantum discord, but spontaneous emission and low detection efficiency have negative impact on the generation of quantum discord. Most importantly, we find that the steady discord depends on feedback parameters, classical driving, and detection efficiency, but not on the initial state.

1. Introduction

In science, the quantification and classification of correlations are of fundamental importance. In the emerging field of quantum information, various correlation metrics have recently been introduced and studied, such as entanglement,[13] quantum correlation, and classical correlation.[4,5] In particular, Ollivier and Zurek[4] and Henderson and Vedral[5,6] have introduced the quantum discord as a measure of quantum correlation beyond entanglement. Recently, quantum discord has been widely explored, including using different alternative methods to guide quantum discord,[7,8] studying quantum discord in different space-time,[9,10] and exploring quantum discord under different interactions.[1114] At the same time, other studies on higher-dimensional and higher spin systems[1518] and experimental advances[19] have also increased interest in quantum discord.

The system considered here consists of two parts, then all correlations can be defined as entropic quantities.[4] In the following, we will define the total correlation by means of measures on one of the subsystems. The quantum system has quantum and classical correlations. All the correlations can be characterized by the quantum mutual information[4,5]

where is the von Neumann entropy. Based on this expression, correlations are generally considered to be separated according to their classical and quantum natures, respectively.[4] In this way, the quantum discord has been introduced as
where is the classical correlation[4,5] defined by the following maximization procedure. A complete set of projector operators must be built for the subsystem B. Then the quantity
must be maximized with respect to the variation of the set of , where , , and .

Since the measurements can be performed on either qutrits, the classical correlation will be calculated by performing measurements on qutrit B. To this end, we consider using the following base vectors for measurement[20]

where (0 ≤ θ, ϕπ/2) and (0 ≤ χ1, χ2 ≤ 2π).

The purpose of this paper is to investigate the effect of a strategic quantum feedback control on the evolution of quantum discord. Quantum feedback is a typical method of controlling quantum system dynamics through external control. Markovian quantum feedback control introduced by Wiseman and Milburn[21,22] has been widely used to feed back measurements to the system to modify the future dynamics of the system, such as suppressing decoherence,[2326] protecting quantum correlations,[27] improving steady state entanglement,[2833] manipulating the geometric phase[34] and stationary state[35] of the dissipative two-level system, preparing a three-atom singlet,[36] W state, and GHZ state,[37,38] and enhancing the accuracy of parameter estimations.[39] Recently, quantum-jump-based feedback control has been used to accelerate quantum evolution,[40] enhance exciton transmission,[41] study on the dissipative stabilization of multipartite entanglement with Rydberg atoms,[42] maneuver Einstein–Podolsky–Rosen Steering,[43] and so on. According to different measurement strategies, homodyne-based and quantum-jump-based feedback controls are proposed. Reference [31] shows that the latter is superior in protecting entangled states. In this paper, we will study the quantum discord of two qutrits under quantum-jump-based feedback control.

In the following, we will study the quantum discord in both feedback and no feedback regimes. The structure of this paper is arranged as follows: In Section 2, we will investigate the evolution of the quantum discord between two qutrits independently coupled to their own environment under quantum-jump-based feedback control. In Section 3, we will investigate the evolution of the quantum discord between two qutrits simultaneously coupled to one environment under quantum-jump-based feedback control. Finally, we give a concise conclusion and discussion.

2. Protection of quantum discord
2.1. Effective master equation

We consider a system consisting of a qutrit which couple resonantly to a damped cavity and is driven by two classical fields. The atomic levels and transitions are depicted in Fig. 1(a). The transition between the levels |3⟩ ← |2⟩ (|2⟩ ← |1⟩) is coupled to the cavity mode a resonantly with the coupling strength g, driven by the classical field with the Rabi frequency Ω. The atom can spontaneously decay from |3⟩ to |2⟩ (|2⟩ to |1⟩) with rate γ1 (γ2), and the cavity mode is damped with rate κa. The complete master equation describing dynamics of our system reads

here, is a superoperator given as
and the superoperator acting on an operator c is defined as . The amplitude damping operator is defined as J = σ12 + σ23, where σkl = |k⟩⟨l|. In what follows, the effective master equation will be derived under the condition of a strongly dissipative cavity, which can be obtained by considering the dynamic evolution of each matrix element of the density operator in the system. In the limit of large decay rate, highly excited modes can be safely ignored, so we only need to expand the density matrix ρ in terms of the a mode number states with small photon numbers.[18,28,36] Therefore, we can get the density operator ρ as
Substituting the above equation into Eq. (4), we can obtain a set of coupled equations of motion for the field-matrix elements

Fig. 1. Schematic view of the model. A qutrit is trapped in a cavity. When a photon from the leaky cavity is detected by a detection apparatus, a feedback pulse is triggered.

Under the condition κag, we make the assumption that and substitute the corresponding results into the expressions of and . After neglecting the terms proportional to second order , we can obtain

Solving the above two equations and adiabatically eliminating[18,28,36,44] the elements ρ2, we can obtain the effective master equation prompting the evolution of the atom
where Γ = 4g2/κa describes the effective decay rate of the atom and the cavity mode a. Moreover, if the collective decay rate Γ is much larger than the atomic spontaneous emission rate γ1(2), we can ignore the effect of γ1(2). The final effective master equation of the system can be expressed as

2.2. Implementation of the feedback

Now let us consider the system with feedback, as shown in Fig. 1(b). The output from the leaky cavity is measured by a detector D, then the measurement is used to trigger a feedback Hamiltonian Hfb = I(t)B. For a feedback scheme based on photodetection measurement, the stochastic master equation for the qutrit system can be described by

where U = exp[−i λ σx] is feedback control with the control strength λ, where . For simplicity, we have assumed γ1 = γ2 = γ.

We consider a system formed by two non-interacting parts A and B. The consistency of each part is shown in Fig. 1. Under this condition, the master equation for the two-qutrit system becomes

In this section, we first assume that two qutrits are in a particular quantum state, then substituting this quantum state into Eq. (12). Finally, we use Eq. (2) to numerically simulate the evolution process of the quantum discord under quantum feedback control.

Now we discuss the case where the initial state is . In Fig. 2(a), we plot the quantum discord of the two qutrits as a function of the dimensionless quantity Γ t and the feedback parameter λ. It can be seen that the evolution of the quantum discord is different under different feedback parameters, and the quantum discord is a periodic function of λ with period T = 2π. Meanwhile, we also find that the decay of the quantum discord can be retarded when some appropriate feedback parameters are chosen. As can be seen from Fig. 2(b), compared with the case of spontaneous emission effects, the quantum discord in the case without spontaneous emission effects has a smaller decay rate. Consequently, the spontaneous emission is a detrimental factor for our current model.

Fig. 2. Quantum discord for the initial state |Ψ⟩ as a function of the dimensionless quantity Γ t and λ with Ω = 0, (a) without (γ/Γ = 0) and (b) with (γ/Γ = 0.01) spontaneous emission effects.

In order to demonstrate the dependency of on classical driving strength Ω, we present Fig. 3, where the initial state is |Ψ⟩ and the plot is drawn at λ = 0.4π. We here choose Ω = 0, Ω/Γ = 2.5, and Ω/Γ = 5. It is easy to find that the classical driving plays a negative role in protecting quantum discord, and as the classical driving strength increases, the impact of the quantum feedback on the dynamics of the quantum discord is complicated. That is to say, our numerical results show that classical driving cannot protect the quantum discord when the quantum-jump-based feedback control exists. In the following, we will focus on the evolution of the quantum discord without classical driving.

Fig. 3. Time evolution of the quantum discord, with the initial state being |Ψ⟩, with λ = 0.4π, for the cases of Ω/Γ = 0, Ω/Γ = 2.5, and Ω/Γ = 5.

Figure 4 illustrates the quantum discord for different initial states (|Ψ⟩, |Φ⟩, and |ϕ⟩) as a function of the dimensionless quantity Γ t, where , , , Ω = 0, and λ = 0.4π. What is shown in Fig. 4 is the quantum discord protection not only depends on the feedback parameter λ and classical driving strength Ω but also on the initial state chosen. At the same time, we use the research methods in Figs. 2 and 3 to study the initial quantum states |Φ⟩ and |ϕ⟩. It is found that the feedback parameters, spontaneous emission, and classical driving have similar effects on the evolution process.

Fig. 4. Evolution of the quantum discord of an initial state |Ψ⟩, |Φ⟩ or |ϕ⟩ with Ω = 0 and λ = 0.4π.

The detection efficiency is another important factor affecting the feedback control process. In order to consider this factor, we can modify the corresponding master equation as

where η represents the detection efficiency, and (1 − η) corresponds to the case where the detector is not clicked and the control action is not applied. In Fig. 5, we plot the quantum discord versus the effect of different detection efficiencies with Ω = 0, λ = 0.4π and two qutrits initially in the state |Ψ⟩. It is shown that the spontaneous emission effect plays a negative role in protecting quantum discord, and the detection efficiency largely determines the protection effect of the quantum-jump-based feedback control on the quantum discord. The protection of the quantum discord is the best when the detection is perfect and there is no spontaneous emission effect. Similarly, we find that other two-qutrit states (|Φ⟩ and |ϕ⟩) show similar behavior to |Ψ⟩ under the appropriate parameters.

Fig. 5. The quantum discord versus the effect of different detection efficiencies with Ω = 0, λ = 0.4π and two qutrits initially in the state |Ψ⟩.
3. Preparation of quantum discord

As can be seen from the above discussion, for the model in Section 2, quantum-jump-based feedback control can protect the quantum discord, that is, reduce the attenuation rate of the quantum discord, but the quantum discord can only be maintained for a limited time. So how do we prepare the quantum discord of two-qutrit system and make it stable? We will discuss it in this section.

3.1. Effective master equation

We consider another system consisting of two qutrits coupled resonantly to a damped cavity. The transition between the levels |2⟩ ↔ |1⟩ and |3⟩ ↔|2⟩ are coupled to cavity mode a resonantly with the coupling strength g, driven by classical fields with the Rabi frequency Ω. The atom can spontaneously decay from |3⟩ to |2⟩ (|2⟩ to |1⟩) with rate γ1 (γ2), and the cavity mode is damped with rate κa.

The complete master equation describing dynamics of our system can be written as

here, is a superoperator given as
and the superoperator acting on an operator c is defined by . The amplitude damping operator is defined as , where σkl = |k⟩⟨ l|. In the limit of large cavity dissipative rates, we can neglect highly excited modes safely.[18,28,36,44] Similar to the analysis procedure of Subsection 2.1, we can obtain the stochastic master equation

3.2. Implementations of the feedback

Now let us consider the system with feedback, as shown in the following Fig. 6.

Fig. 6. Schematic view of the model. Two qutrits are trapped in a cavity. When a photon from the leaky cavity is detected by a detection apparatus, a feedback pulse is triggered.

The output from the leaky cavity is measured by a detector D, then the measurement is used to trigger a feedback Hamiltonian Hfb = I(t)B. For a feedback scheme based on photodetection measurement, the stochastic master equation for the two-qutrit system can be written as

where , with . For simplicity, we have assumed . We choose the local feedback that, like the spontaneous emission term, cannot maintain the symmetry of the atomic exchange. This can allow feedback to move the population between subspaces, limiting the likelihood of destructive interference.

The quantum discord as a function of feedback and driving strength is shown in Fig. 7. Compared with the case of spontaneous emission effects, the quantum discord of the composed system in the case without spontaneous emission effects has a large maximum value in steady state, therefore we hold the point that the atomic spontaneous emission is a detrimental factor for our current model. And it is obvious that a nearly perfect quantum discord can be achieved by quantum-jump-based feedback.

Fig. 7. The quantum discord as a function of feedback and driving strength (a) without (γ/Γ = 0) and (b) with (γ/Γ = 0.01) spontaneous emission effects.

In order to demonstrate the dependency of the quantum discord on classical driving strengths Ω, we present Fig. 8, where the initial state is |Ψ⟩ and the plot is drawn at λ = π. We here choose Ω = 0, Ω/Γ = 2.5, and Ω/Γ = 5. We can see from Fig. 8 that when the feedback parameter λ = π is chosen, the effect of classical driving on the dynamics of the quantum discord is complicated. But it is easy to find that the classical driving plays a positive role in the preparation of the quantum discord, and with the advent of classic driving, a nearly perfect maximum value of discord appears. That is to say, our numerical results show that when the appropriate quantum-jump-based feedback control is present, the classical driving is very beneficial for the preparation of the quantum discord.

Fig. 8. Time evolution of the quantum discord, with the initial state being |Ψ⟩, with λ = π, for the cases of Ω/Γ = 0, Ω/Γ = 2.5, and Ω/Γ = 5.

Figure 9 illustrates the quantum discord for different initial states (|Ψ⟩ and |Φ⟩) as a function of the dimensionless quantity Γ t, where , |Φ⟩ = |11⟩, Ω/Γ = 1, and λ = π. Figure 9 shows that the quantum discord evolution process not only depends on the feedback parameter λ and the classical driving strength Ω but also on the initial state chosen. For any initial state, the long-term evolution of the quantum discord can reach the maximum.

Fig. 9. Evolution of the quantum discord of an initial state |Ψ⟩ or |Φ⟩ with Ω = Γ and λ = π.

The detection efficiency is another important factor affecting the feedback control process. In order to consider this factor, we can modify the corresponding master equation as

where η represents the detection efficiency, and (1 − η) corresponds to the case where the detector is not clicked and the control action is not applied. In Fig. 10, we plot the quantum discord versus the effect of different detection efficiencies with Ω = Γ, λ = π, and two qutrits initially in the state |Ψ⟩. Obviously, the system always achieves the maximum value of discord without atomic decay, but the time required depend on the efficiency of the detection process. Furthermore, the presence of spontaneous emission leads to a reduction in the final quantum discord, which is also affected via the inefficiency of the detector, this conclusion is consistent with the conclusion in the literature.[36]

Fig. 10. The quantum discord versus the effect of different detection efficiencies with Ω = Γ, λ = π, and two qutrits initially in the state |Ψ⟩. Without decay the system always reaches the maximum value of discord, with the time depending on the efficiency, while the presence of spontaneous emission results in a decrease in the final quantum discord.

To consider the correctness of the adiabatic approximation, we can rewrite the main equation (14) as

In the adiabatic regime with κa = 1600γ, g = 200γ, Ω = 100γ, and λ = π, the discord of the steady state is about 1.39 that is close to the result of effective model denoted by the red line in Fig. 10.

4. Discussion and conclusions

In general, the existence of decoherence makes it an obstacle to the quantum information science to experimentally produce faithful and reliable entangled states, because quantum system will inevitably interact with their surrounding and cause information loss. In the context of cavity quantum electrodynamics, the typical decoherence factors include spontaneous emission of atoms and cavity decay. In the absence of atomic spontaneous emission, the reason why quantum feedback can protect quantum discord is that the fast damping rate κa of the cavity makes it easy for photon to leak from the cavity rather than interact with the atom. Once the detector detects a leaked photon, the feedback control will act on the system to modify the future dynamics of the system. In this way, proper feedback control suppresses the decoherence of the system and the surrounding environment, and protects the quantum discord of the quantum system. The presence of spontaneous emission will aggravate the loss of information in the quantum system and have a negative impact on the protection of quantum discord. In the preparation of quantum discord, we adopt the basic idea of direct quantum feedback, that is, in a system consisting of largely decay cavity and negligible atomic spontaneous emission, the feedback with not preserving the symmetry with respect to exchange of atoms can be used to design an approximate maximum discord steady state. When the spontaneous emission of an atom occurs, the steady state quantum discord value decreases. It is due to a competition between the spontaneous emission and the feedback rate.[18,29] So the atomic decay still plays a negative role in the current scheme.

In experiments, for an optical cavity approaching the bad cavity limit, where the atom–cavity cooperativity parameter is C = g2/κa γ ≈ 19[45] with (κa, g, γ) ∼ 2π × (255.05, 70, 1) MHz,[46] the quantum discord is about 1.24 based on Eq. (19) in the case of Ω = g and λ = π. The value is relatively high for the discord. In the future cavity quantum electrodynamics experiments, it is expected to obtain a better bad cavity structure. The atomic configuration involved in our scheme may be implemented with metastable Krypton, whose relevant atomic levels are shown in Ref. [47]. However, its structure does not fully satisfy our theoretical assumptions. Therefore, we also hope that in the near future, the cascaded three-level atom that fully satisfies our hypothesis can be found.

In summary, we find that a quantum-jump-based feedback scheme can protect and prepare the quantum discord of qutrit-qutrit system. In Section 2, we show that quantum feedback can protect the quantum discord, and the emergence of spontaneous emission and classical driving are negative factors. In addition, we find that the feedback has different protection effects on the quantum discord of different initial states, and the protection effect is the best for the quantum discord in the initial state |Ψ⟩. Finally, we discuss the effects of low efficiency and spontaneous emission on quantum discord during the detection process. We find that high detection efficiency is more conducive to the protection of the quantum discord, and the existence of spontaneous emission will undermine the protection of the quantum discord. In Section 3, we show that another quantum-jump-based feedback scheme can produce stable maximum quantum discord. First, the quantum jump-based feedback scheme can improve the steady-state discord of the two-qutrit system, but the spontaneous emission in the current scheme plays a negative role. In addition, we explore the effects of classical driving on quantum discord and find that the existence of classical driving is a positive factor and a necessary condition for preparing the maximum steady-state discord. Later, we find that the feedback can produce a stable maximum quantum discord, regardless of the initial state. Finally, we discuss in detail the effects of low efficiency and spontaneous emission on quantum discord during the detection process. On the one hand, the system always achieves the maximum value of discord without atomic decay, but the time required depends on the efficiency of the detection process. On the other hand, the presence of spontaneous emission leads to a reduction in the final quantum discord, which is also affected by the inefficiency of the detector. Therefore, in the two-qutrit system, quantum-jump-based feedback is an effective scheme for protecting and preparing quantum discord.

Reference
[1] Bennett C H DiVincenzo D P Smolin J A Wootters W K 1996 Phys. Rev. 54 3824
[2] Vidal G 2000 J. Mod. Opt. 47 355
[3] Horodecki R Horodecki P Horodecki M Horodecki K 2009 Rev. Mod. Phys. 81 865
[4] Ollivier H Zurek W H 2001 Phys. Rev. Lett. 88 017901
[5] Henderson L Vedral V 2001 J. Phys. A: Math. Gen. 34 6899
[6] Vedral V 2003 Phys. Rev. Lett. 90 050401
[7] Shi W N Wang D Sun W Y Ming F Huang A J Ye L 2018 Laser Phys. Lett. 15 075202
[8] Chen M N Sun W Y Huang A J Ming F Wang D Ye L 2018 Laser Phys. Lett. 15 015207
[9] Ming F Wang D Liu Y 2019 Ann. Phys. 531 1900014
[10] Wang D Shi W N Hoehn R D Ming F Sun W Y Kais S Ye L 2018 Ann. Phys. 530 1800080
[11] Li H Z Han R S Zhang Y Q Chen L 2018 Chin. Phys. 27 120304
[12] Faizi E Eftekhari H 2015 Chin. Phys. Lett. 32 100303
[13] Zheng Y D Mao Z Zhou B 2018 Chin. Phys. 27 090306
[14] Zad H A 2016 Chin. Phys. Lett. 33 090302
[15] Guo K T Xiang S H Xu H Y Li X H 2014 Quantum Inf. Process. 13 1511
[16] Maleki Y 2016 Quantum Inf. Process. 15 4537
[17] Lopez C E Lastra F 2017 Phys. Rev. 96 062112
[18] Wu Q C Ji X 2013 Quantum Inf. Process. 12 3167
[19] Dakic B Lipp Y O Ma X S Ringbauer M Kropatschek S Barz S Paterek T Vedral V Zeilinger A Brukner Č Walther P 2012 Nat. Phys. 8 666
[20] Caves C M Milburn G J 2000 Opt. Commun. 179 439
[21] Wiseman H M Milburn G J 1993 Phys. Rev. Lett. 70 548
[22] Wiseman H M 1994 Phys. Rev. 49 2133
[23] Viola L Lloyd S 1998 Phys. Rev. 58 2733
[24] Katz G Ratner M A Kosloff R 2007 Phys. Rev. Lett. 98 203006
[25] Ganesan N Tarn T J 2007 Phys. Rev. 75 032323
[26] Zhang J Li C Wu R Tarn T Liu X 2005 J. Phys. A: Math. Gen. 38 6587
[27] Yu M Fang M F 2017 Int. J. Theor. Phys. 56 1937
[28] Wang J Wiseman H M Milburn G J 2005 Phys. Rev. 71 042309
[29] Carvalho A R R Hope J J 2007 Phys. Rev. 76 010301(R)
[30] Li J G Zou J Shao B Cai J F 2008 Phys. Rev. 77 012339
[31] Carvalho A R R Reid A J S Hope J J 2008 Phys. Rev. 78 012334
[32] Chen Y Zou J Li J G Shao B 2010 Acta Phys. Sin. 59 8365 in Chinese
[33] Chen L Wang H F Zhang S 2014 Chin. Phys. 23 030301
[34] Sun H Y Shu P L Li C Yi X X 2009 Phys. Rev. 79 022119
[35] Wang L C Huang X L Yi X X 2008 Phys. Rev. 78 052112
[36] Shao X Q Zheng T Y Zhang S 2012 Phys. Rev. 85 042308
[37] Sun W M Su S L Jin Z Liang Y Zhu A D Wang H F Zhang S 2015 J. Opt. Soc. Am. 32 1873
[38] Chen L Wang H F Zhang S 2013 J. Opt. Soc. Am. 30 475
[39] Zheng Q Ge L Yao Y Zhi Q J 2015 Phys. Rev. 91 033805
[40] Yu M Fang M F Zou H M 2018 Chin. Phys. 27 010303
[41] Ji Y Q Qin M Shao X Q Yi X X 2017 Phys. Rev. 96 043815
[42] Shao X Q Wu J H Yi X X 2017 Phys. Rev. 95 022317
[43] Huang Z M 2018 Int. J. Theor. Phys. 57 3473
[44] Wiseman H M Milburn G J 1993 Phys. Rev. 47 642
[45] Metz J Beige A 2007 Phys. Rev. 76 022331
[46] Dayan B Parkins A S Aoki T Ostby E P Vahala K J Kimble H J 2008 Science 319 1062
[47] Pu H Cai T Bigelow N P Grove T T Gould P L 1995 Opt. Commun. 118 261